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We review the entanglement properties in collective models and their relationship with quantum 
phase transitions. Focusing on the concurrence, which characterizes the two-spin entanglement, 
we show that for first-order transitions, this quantity is singular but continuous at the transition 
point, contrary to the common belief. We also propose a conjecture for the concurrence of arbitrary 
symmetric states which connects it with a recently proposed criterion for bipartite entanglement. 
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I. INTRODUCTION 

The relationship between entanglement properties and 
quantum phase transitions (QPTs) has been recently the 
subject of many debates. An important issue is to under- 
stand how the structure of the ground state is affected by 
the existence of a critical point. The first answers came 
with the study of one-dimensional spi n sy stems QJ 
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where a strong interplay between QPTs and the variation 
of the entanglement with the control parameters has been 
in evidence. However, a complete classification of the 
possible scenarios likely to arise depending on the nature 
of the transitions is still lacking [22, . In this perspec- 
tive, collective models provide a good framework to inves- 
tigate these questions. Indeed, as shown in the Lipkin- 
Meshkov-Ghck modeljpllllllllllll, the Dicke model 
of superradiance ^,"3% or the two-level BCS model 
the essential features met in low-dimensional systems are 
present in these systems which can be viewed as infinite 
dimensional. The main advantage of these models is that 
their phase diagram can be determined within a mean- 
field approach and that their entanglement properties, 
encoded in finite-size corrections, can be computed with 
standard methods. 

Characterizing the entanglement of a quantum state is 
not an easy task and there actually exists only a few re- 
liable tools. For many-body states, one only knows how 
to quantify bipartite entanglement through various mea- 
sures such as the block entropy (obtained after a partition 
of the system into two arbitrary blocks) or the concur- 
rence l32| which characterizes the entanglement between 
two spins after tracing out over all others. For collec- 
tive models, this latter quantity is particularly well-suited 
since it does not depend on the two spins selected, all 
spins being completely equivalent. 

The goal of the present study is two fold. First, 
we discuss the concurrence for symmetric (permutation- 
invariant) states. The motivation for considering these 
states is that the ground states of collective models are 
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(often) invariant under the permutation group. Simple 
expressions of the concurrence have already been ob- 
tained for some particular such states (parity symmetric) 
{33i | but we propose here an extension of these results for 
a wider class. This extension leads us to conjecture a 
general form of the concurrence for arbitrary symmet- 
ric states that we relate to a recently proposed criterion 
[33. 135I I . Second, we analyze the behavior of the concur- 
rence in several models in which first- and second-order 
transitions occur. The main result of this paper is that 
there exist situations where a discontinuity of the order 
parameter (associated with a first-order QPT) does not 
imply a jump of the concurrence at the transition point. 
This is a rather surprising result since the order parame- 
ter appears explicitly in the reduced density matrix from 
which the concurrence is extracted. 

This paper is organized as follows. In Sec. ^ we de- 
rive a simple formula for the concurrence of a class of 
symmetric states and propose a conjecture for arbitrary 
symmetric states. In Sec. IIIII we introduce the collec- 
tive models considered in the next sections. Section HVI 
is devoted to a review of the concurrence in collective 
models in perpendicular magnetic field which are known 
to display second-order QPT. In Sec. we consider the 
most general uniaxial collective model. After discussing 
its phase diagram, we compute exactly the concurrence of 
the ground state in the thermodynamical limit as a func- 
tion of the magnetic field and show that it is always con- 
tinuous at the transition point. Finally, we conclude and 
give some perspectives in Sec. IVII In the Appendix, we 
establish the correspondence between the uniaxial model 
in an arbitrary field and a two-level boson model recently 
proposed in nuclear physics which allows a simple 
physical interpretation of the order parameter. 



II. A CONJECTURE FOR THE 
CONCURRENCE OF SYMMETRIC STATES 

In this section, we give some simple expressions for the 
concurrence for symmetric states and conjecture a rela- 
tionship with a recently derived bipartite entanglement 
criterion 

The concurrence yJl is defined as follows. Let us con- 
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sider a system of N spins half described by a density ma- 
trix. Let be the reduced density matrix associated 
with spins i and j obtained by tracing out the density 
matrix over aU spins except spins i and j. Next, let us 
introduce the spin-flipped density matrix 



Pr,3 = i'^y ^ '^y) Ph (o-y ^^y), 



(1) 



where p* j is the complex conjugate of pij. The concur- 
rence related to spins i and j is then defined by 



Cij — max |o, \/Xi — y/X2 — 



A,- 



A4}, (2) 



where the A^'s are the four real positive eigenvalues of 
PijPij, and where Ai > Ai+i. If Cij > 0, spins i and j 
are said to be entangled. 

For symmetric states, it is clear that pij and hence Cij 
do not depend on the specific choice of i and j so that we 
will omit these indices in the following. Now, let us try 
to find a simple, and physically meaningful, expression 
of the concurrence. As shown in Ref. 33], the reduced 
density matrix for symmetric states in the standard basis 
{I TT),| n),| it), I ii)} (with a,| T,i) =±l T,i)) can 
generically be written as: 



X 



+ 



X, 



with 



v± 



x± 



u 

x+ w y x*_ 
x+ y w xl 

U X- X- V- 



N^-2N + 4(5f) ± 4{N - 1)(S,) 
4:N{N - 1) ' 
iN~l){S+)±{[S+,S.U) 
2N{N-1) 

4:N{N - 1) ' 



(3) 



iV(iV- 1) 
I) 



7V(iV- 1)' 



(4) 
(5) 
(6) 
(7) 
(8) 



where [A, = AB+BA. Note that since for symmetric 
states = ■j(y + 1) one has w = y. 

The only simple expressions of the concurrence pro- 
posed in the literature concern symmetric states that 
are spin-flip (or parity) symmetric, i.e., that are eigen- 
states of the operator 11 = Ilicr*. For such states, one 
has x± — 0, and the concurrence is given by 33] 



2max{0, \u\ — y} 



if 2y < y/v+v- + \u\ 



2 max{0, y — y/v+v-} if 2y > y/v+v- + \u\ 



(9) 

An important remark is that, for real states, this result 
can also be written as 



C = TnaK{0,Cy,Cz}, 



(10) 



where, for a given direction n, we have introduced the 
quantity 



r 

^71 



1 



2N{N - 1) 



(11) 



N{N-2) + A(S?^) ~ 4{N-l)(Sn) 



Note that for eigenstates of H, one has (Sy) — so 
that for n^y, one has {N ~ l)Cy = 1 - 4(3^) /N. 

Before coming to our conjecture, let us also compute 
the concurrence for another type of states which break 
the parity (spin-flip) symmetry but which have real co- 
efficients in the standard basis {\N/2, AI)} of and Sz- 
For such states, one has 



{Sy) = {[S.,Sy] + ) = {[Sz,Sy] + )=Q, 



(12) 



SO that p ^ p* and the characteristic polynomial associ- 
ated with pp factorizes: 

P(A) = det(pp-A), (13) 
= det (ptTy ® ay — a/a) det [pay ® ay + V\) , (14) 
= Q{p)Q{~p), (15) 

where we have set p = V\. The A.^'s of pp are thus given 
by the squares of the roots of Q. Further, since w = y, 
one has P(0) = 0. We thus have to analyze the roots of a 
third-order polynomial. Denoting by pi the three other 
roots of Q, one has 



Q{p) 



p 



(16) 



p^ + 2{u- y)p^ + 

p{v? — v^v- — Auy + Ax^x^) + (17) 
2(yu-|_i;_ + 2ux^x^ — v?y — v^x?_ — w_a;^). 



The concurrence can then, under some conditions, be 
directly related to the piS. Indeed, if the piS are such 
that /ii < and < p2,?, < then one has 



C = max |o, — 

i 

= max{0,2(u - y)} , 
= max{0,Cy}. 



(18) 

(19) 
(20) 



The fact that the direction y plays here a special role is 
due to the fact that the states considered have real coef- 
ficients in the basis {\N/2, M)} or, equivalently, satisfy 
(|12|l . Finally, a close analysis of Q allows one to show 
that the conditions pi < and < /i2,3 < \pi\ are sat- 
isfied as soon as u > w. When this last condition is 
violated, we have no explicit form of the concurrence. 

This analysis leads us to following conjecture: the con- 
currence of any arbitrary symmetric state is given by 



; |o, maxC„| 



(21) 
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where Cn is given by Eq. We have checked numer- 

ically by choosing symmetric states with random com- 
plex coefficients in the basis that, indeed, one always has 
max„ Cn — VA]" — ^/\2 — a/As — VXj • Although we have 
no proof of this conjecture, we would like to give a strong 
argument in its favor. Recently, a necessary and sufficient 
criterion for bipartite entanglement in symmetric states 
has been proposed by Korbicz et al. 35]. This cri- 
terion indicates that bipartite entanglement is present if 
and only if there exists a direction n such that 



1 



7V2 



-4 



N 



> 0. 



(22) 



which is completely equivalent to Cn > 0. In this case, 
following Eq. (|^ . C„ > implies C > 0, which is 
indeed a necessary and sufficient condition for bipartite 
entanglement. 

Now, we shall use this simple expression (|21ll of the 
concurrence to characterize the entanglement of the 
ground state in collective models. 



III. GENERALITIES ABOUT COLLECTIVE 
MODELS 

The collective models considered here are defined as 
systems in which all spins mutually interact. The Hamil- 
tonian H of such systems can thus be expressed in terms 
of the total spin operators Sa — J2i '^a/'^ where the ct^'s 
are the Pauli matrices. As a direct consequence, one has 
[H, S^] = 0. If one considers only two-spin interactions 
and a coupling to a magnetic field h, the most general 
collective model Hamiltonian can be written as 



Hn = 



^^lijSiSj + h.S, 



(23) 



The interaction matrix can be diagonalized so that, af- 
ter a rotation and up to some simple redefinitions, this 
Hamiltonian also reads 



Hi 



(24) 



In a fixed spin sector, one further has = 5^ + 5^ + Sl 
so that, finally, the most general Hamiltonian reads 



(25) 



The most general parameter space to be investigated is 
thus, after normalization, four dimensional. In this work, 
we restrict our discussion to ferromagnetic interactions 
lx,y < for which the ground state of H belongs to the 
symmetric space — y(t ~^ 

For such models, a variational approach assuming a 
completely separable wave function is very efficient to 
obtain the ground state in the thermodynamical limit. 
However, if one is interested in the entanglement proper- 
ties, a separable state does not bring much information. 
As discussed in Ref. 01 , if the concurrence of the ground 



state C indeed vanishes in the large- iV limit, the rescaled 
concurrence Cr = (N — 1)C has interesting properties in 
this limit. In the following, we will consider this latter 
quantity and focus on two categories of systems: the bi- 
axial model in a transverse field and the uniaxial model 
in an arbitrary field. 



IV. BIAXIAL MODEL IN TRANSVERSE FIELD 

Let us first consider the simplest nontrivial situation 
where the magnetic field is perpendicular to the inter- 
action directions. Generically, the Hamiltonian of this 
system can be written as 



hzSz 



(26) 



where 7 is the anisotropy parameter. Here, since we re- 
strict our discussion to ferromagnetic coupling, we sup- 
pose, without loss of generality, that < 7 < 1. We 
also suppose in the following that hz > 0. As widely dis- 
cussed in the literature (see, e.g., Ref. I23|); this system 
displays a second-order QPT at /iz = 1. A good order 
parameter for this transition is 



m^l-4{S^,)/N^ 



l-hl 




for 
for 



0<hz<l, 
hz > 1. 



(27) 



For a discussion of the critical exponents and the finite- 
size corrections associated with this transition, we refer 
the interested reader to Refs. [l^li^. In the following, 
we simply recall the main results concerning the concur- 
rence in these models. 



A. The anisotropic case 

As shown in Ref. [23|, the rescaled concurrence for 
< 7 < 1 is given by 



N 



iN-l)Cy, 

(N - l)Cz, 



(28) 
(29) 



where C„ is defined by Eq. Hll|) . These expressions H28|) 
and (|29|l underline the role played by the special point 

hz = ■^^^v~^'\/^ which the two ground states can be 
chosen as separable 

m, 0)) = [cos (0/2) e-'^^/2| t)z+sin (9/2) e^^/^H), 

(30) 

with 9 — arccos hz and = or tt. Such a special point is 
often present in anisotropic models |38j | . Apart from this 
point where C = for all N, the rescaled concurrence 
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can be computed in the thermodynamical limit 



UNIAXIAL MODEL IN ARBITRARY FIELD 



C^^^^ = 1 



1-7^ 



(31) 
(32) 
(33) 



Obviously, when hz goes to infinity, the rescaled con- 
currence vanishes since the ground state is, in this limit 
\ijj{0,0)) (fully polarized state in the z direction). In the 
zero- field limit, the rescaled concurrence is nontrivial ex- 
cept for 7 = where the ground states \tlj{n/2, (j) = 0, tt)) 
(fully polarized states in the x direction) have a vanishing 
rescaled concurrence. 



B. The isotropic case 

In the isotropic case 7 = 1, the problem becomes sim- 
ple since one further has [H, Sz] — and the eigenstates 
of the Hamiltonian arc thus the eigenstates Af)} of 
and Sz- The ground state is obtained for S — N/2 
and the following values of M 



-I{hzN/2) 
-N/2 



for 0<hz <l, 
for hz > 1, 



(34) 



where I{x) gives the integer part of x, in the following 
sense : ii x = X + Sx with X an integer and Sx G [0, 1[, 
then I{x) = X for 5x G [0, 1/2 [ and I{x) = X + 1 for 

fee [1/2, 1[. 

As shown in Ref . [23 j the rescaled concurrence of the 
state \N/2, M) is: 



Cr = {N- \)Cz 



(35) 



Note that these states, known as Dicke states, are par- 
ity symmetric and this expression simply comes from Eq. 
(|10|l . In the thermodynamical limit, one thus has a dis- 
continuity at the critical point hz — 1 where the rescaled 
concurrence jumps from for Mq = —N/2 {hz > 1) to 2 
for Mq = —N/2 -\- 1, and goes to 1 in the zero-field limit, 
i.e., when Mq goes to 0. 

These results show that the rescaled concurrence is sen- 
sitive to the existence of a QPT since it is maximum, and 
singular, at the critical point ft,^ = 1. Further, it distin- 
guishes between the two universality classes 7 7^ 1 and 
7 = 1. Indeed, in the former case, Cr is continuous at the 
transition point whereas in the latter, it displays a jump 
despite the continuous character of the transition. To our 
knowledge, it is a unique example of a second-order QPT 
associated with a discontinuous rescaled concurrence. 



We now turn to the main contribution of this study 
which concerns the uniaxial model in an arbitrary field 
whose Hamiltonian can generically be written as 



H 



-L:|| 

X 



hxSx + hzSz 



(36) 



with hz > 0. As discussed in the last section for h^ = 0, 
this model displays a second-order transition point when 
varying hz with a critical point at /i^ = 1. As we shall 
see, for hz < 1, one also faces a QPT when varying the 
parallel field h^. This transition is first order and oc- 
curs at hx — 0. Within the parametrization H36|l . it is 
not easy to build an order parameter characterizing this 
transition. However, in the Appendix, we give the cor- 
respondence between this model and a two-level boson 
problem recently introduced in nuclear physics which 
allows one to get an order parameter. For /i^ > 1, a sim- 
ple mean-field analysis similar to those presented in Ref. 
|27l | also allows one to show that there is no transition 
when varying h^- Similarly, for /la, 7^ 0, no transition is 
found when varying hz- 

Let us now discuss the behavior of the concurrence in 
the plane {hx,hz). The ground state of ^f^'" is in the 
maximum spin sector and has real coefficients in the basis 
{|A^/2,M)}. Further, as we shall show, one has Cy > 
so that, as discussed in Sec. ^ its rescaled concurrence 
is given by 



Cr = (iV - l)Cy 



(37) 



To compute this quantity in the thermodynamical 
limit, we follow the same line as for the transverse field 
case ■27^. The first step is to introduce the Holstein- 
Primakoff representation of the spin operators: 



a - N/2, 



S+ = iVi/V(l-7iJiV)i/' = (5_)t, 



(38) 
(39) 



with S± = Sx ±iSy and Ua = a' a. The creation (a' ) and 
annihilation (a) bosonic operators satisfy the canonical 
commutation rules [a,a^] = 1. The standard way to 
expand these operators is to suppose {na)/N <C 1, i.e., 
to assume that the ground state, in the thermodynamical 
limit, is fully polarized in the z direction. In the opposite 
case, one can either perform a rotation to bring the z- 
axis along the semiclassical spin direction, or shift the 
bosonic operators by setting a} = y/N A -I- 6^ . This latter 
procedure provides a macroscopic expectation value of 
Sz which is of order N, and one then has {ni,)/N ^ 1. 
We adopt this latter approach but, of course^the former 
is strictly equivalent. As explained in Ref. j27l |. we need 
to expand the Hamiltonian only at order (1/A^)'^ to get 
the rescaled concurrence in the thermodynamical limit. 
At this order, the Hamiltonian ()36|l reads 



5 



H 



X 



2(l + /32)2 



2(1+^2)3/2 



1 



^ I V" 8/3^ + /i,/3(2 + 3/32 + /34)] + 2nfc [2(1 - 6/3^) + /i,/3(4 + 5/3^ + /34) - 4/i,(l + /S^)] 

(40) 



+2 + /32 [ - 4 + KI3{1 + /32)] } + 0{\/^/N), 



where, for convenience, we have set A = P/y/TVW- To 
diagonahze i?3^'" at this order, we proceed in two steps. 




FIG. 1; (Color online) /3o as a function of hx and h^. The 
singular behavior is clearly observed &t = G for < 1. 
For > 1, /3o is a smooth function for any h^- 

First, we choose /3 such that it minimizes the classical 
energy E{f3, h^, h^), i.e., the term proportional to N. 
This choice, /3o, cancels the term of order ^/N. Then, 
one simply has to diagonalizc a quadratic form, which is 
straightforward. A simple study of the classical energy 
shows that (3o is an odd function of for all 7^ 0. For 
= 0, one has /3§ = (1 — hz)/{l + h^) for < I and 
(3q = otherwise. As a by-product, contrary to /3o, is 
a continuous function of the parameters and h^. As 



explained in the beginning of this section, several regimes 
must thus be distinguished. 

1. ft.^ < 1: /3o is discontinuous at ~ signaling 
the presence of a first-order QPT on the line {h^ — 
0,h, < 1) ; 

2. hz — I: /3o is continuous at — but dp^/dhx 
and df3o/dhz are discontinuous functions at the 
point {hx = 0,hz — 1). One thus has a second- 
order transition point there ; 

3. /iz > 1: /3o is a regular function and so that 
no transition is expected in this region. 



For illustration, we plotted /3q as a function of and hz 
in Fig. H] 

Once /3o is determined we diagonalize the quadratic 
Hamiltonian H43|l via a standard Bogoliubov transform 
by setting 

= cosh(e/2)6^ + sinh(e/2)6 = (c)^ (41) 

where O is such that the Hamiltonian expressed in terms 
of the c's is diagonal. This leads us to choose 

, ^ _ -2 + 8/3g - KPo (2 + 3/3g + Pl) 

' 2(1 - 6/32) + h^PoiA + 5/32 + /34) - AKil + Pi) ■ 

(42) 

Finally, one obtains 



rrUl _ ,, -2/30^ + 2fe./3o(l + P',) - hz{l - P^) 

^ " 2{l + Pir 

2^ I fl„('_9fln -^ h-.n + /32^ "1/2 



where 



^ 2-12/3o2-4/^, l + /3o' +/i,/3o 4 + 5/3o'+/3o^ - 2 - 8/3o' + /i,/3o 2 + 3/3o' + /3o^ ^^^^ 

=(/3o, /iz) = ^ ^ ^ —■ (44) 

16(1 + /3o') 



r 



The ground state, at this order, is given by the zero c-boson state so that we now need to write the rescaled 
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concurrence for the ground state in terms of these bosons. 
Using Eqs. (EH), igTl), and lO, one finds 

Or — i ^ . 

1 + /3o' 

(45) 

For hx — 0, one retrieves expressions ()31|l and (|32|l set- 
ting 7 = 0. The most interesting result concerns the 
behavior of Cr at the transition point — 0. Indeed, 
as explained above, /3q is a continuous function of and 
hz so that Cr is also a continuous function of these pa- 
rameters. In particular, when hz < I, the system under- 
goes a first-order QPT but the concurrence is continuous 
even though the order parameter has a jump. However, 
it can easily be checked that for = and < hz < 1, 
dC^i/dhx is discontinuous but finite whereas dC^/dhz is 
continuous (see Fig. By contrast, for h^ — and 
hz = I, dCn/dhx and dCn/dhz are discontinuous and 
diverge. Finally, for /i^ > 1, Cr is a smooth function of 
hx and hz as can be seen in Fig. El 




FIG. 2: (Color online) Rescaled concurrence as a function of 
and hz in the thermodynamical limit. Despite the presence 
of a first-order transition line (for h^ = G and /iz < 1), Cr is 
a continuous function of these parameters. 

This example underlines the fact that the behavior of 
this entanglement measure does not directly inform us 
about the nature of the transition. It is likely that a 
complete classification certainly requires more informa- 
tion as recently discussed in Ref. [23. This leads us to 
comment on the theorem proposed in Ref. p3 | which 
states that a discontinuous concurrence is necessary, and 
sufficient, to signal a first-order QPT. Our result does 
not violate this theorem since one of its assumptions is 
not satisfied. Indeed, we face here a situation where the 
discontinuous elements of the reduced density matrix p 
(for example {Sx)) do not appear in the concurrence. 

Another interesting point that has attracted much at- 
tention recently concerns the finite-size corrections to the 
concurrence. As already discussed in Refs. [2^ |2^ 
for hx = 0, the only divergence in the 1/A'' expansion of 
Cr occurs for /i^ = 1, detecting a finite-size scaling expo- 
nent for Cr which is 1/3. For hx = and hz < 1, we have 



checked that the 1/A'' correction to {Sy)/N is regular so 
that one expects Cr to have simple 1/N corrections in 
the large- A'' limit. 



VI. CONCLUSION AND PERSPECTIVES 

In this work, we proposed a conjecture for the concur- 
rence of arbitrary symmetric states. This conjecture has 
been checked numerically and is deeply connected with 
a recently proposed criterion for bipartite entanglement 
|35l |. Using this conjecture, we have analyzed the en- 
tanglement properties of a collective model with uniaxial 
interaction and arbitrary field which displays first- and 
second-order QPTs. For this system, we have exactly 
computed the rescaled concurrence Cr in the thermo- 
dynamical limit and in the whole parameter range. Al- 
though Cr is singular at the transition point, it is never 
discontinuous at the transition point as one could have 
expected. 

The next step would be to consider the most gen- 
eral (ferromagnetic) collective system which is the biaxial 
model in arbitrary field given by the Hamiltonian 

H^-j = -^{Sl+jS^y) + h.S. (46) 

In this case, the parameter space to investigate is four di- 
mensional and the eigenstates have complex coefficients 
in the standard basis {|A^/2,M)} so that we cannot, a 
priori, analytically determine the direction n giving the 
concurrence via our conjecture. Such a (numerical) study 
is beyond the scope of the present work but would cer- 
tainly be of interest to understand the competition be- 
tween the interaction anisotropy and the orientation of 
the magnetic field. 

It would also be instructive to study the one- 
dimensional XY model in arbitrary field which is the 
counterpart of (|46|l in finite dimensions. Even in the uni- 
axial (Ising) case, the behavior of the concurrence as a 
function of the field is an interesting issue and a com- 
parison with the results of the present work would bring 
complementary informations about the interplay between 
entanglement and QPT. 

Finally, one may also think about using other entan- 
glement measures such as the block entropy. For collec- 
tive models in perpendicular field, it has already revealed 
nontrivial behavior l^^. A complete study of the en- 
tanglement entropy for the ground state of H^^y, would 
definitely be enriching. 
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APPENDIX A: CORRESPONDENCE BETWEEN 
A TWO-LEVEL-BOSON MODEL AND THE 
UNIAXIAL MODEL IN ARBITRARY FIELD 



with tana = — y/2. Such a rotation diagonalizes the 
interaction matrices so that the Hamihonian (|A1(I is now 
given by: 



In a recent paper [3g , a simple two-level boson model 
has been introduced. The Hamiltonian of this model is 



H 



1 



X nt 



N 



(Al) 



where the operators rit and are defined as 

nt = th, Qy = sh + t^s +y th, (A2) 

in terms of two species of scalar bosons s and t, x and 
y being two independent control parameters. The total 
number of bosons iV = + nt is a conserved quantity. 
To map the Hamiltonian (|A1|) onto the Hamiltonian of 
the uniaxial model, we use the Schwinger representation 
of the spin operators: 



(A3) 



In terms of the spin operators, the Hamiltonian ljAl|l then 
reads 



H = x[S 



1 



N 



25, 



+ y(5. + f)]'. (A4) 



Next, let us perform a rotation around the y axis: 



COS a sin a 

1 
— sin a cos a 




(A5) 



N 



H = X Ny 

2 ^ A 



1-x 52 (l"^)(4 + 2/')' 



AN 



cos a 



S^xcosa + S j^^^y^l^'' ""'K osa, (A6) 



where Xc = (4+y^)/(5 + 2/^) is known to be the transition 
point [3^. The correspondence with the uniaxial model 
Hamiltonian H36|l is then straighforward. Indeed, after 
normalizing the Hamiltonian (|A6() to have an interaction 
term which is —S'^/N the mapping reads 



y{5 + y'^){x - Xc) 
(a; - 1) (4 + 2/2)3/2' 

2x 

(2;- 1) (4 + 2/2)3/2- 



(A7) 
(A8) 



One can check that the transition point x — Xc corre- 
sponds, as expected, to hx = 0. Let us also mention that 
in order to keep a ferromagnetic interaction term, one 
must have x < 1 which implies some constraints on the 
field amplitude. 

The main advantage of the parametrization ljAl|l . 
in terms of the two bosons s and t, is that the order 
parameter is simply (nt)/N j36j |. 
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